Coupled responses of thermomechanical waves in functionally graded viscoelastic nanobeams via thermoelastic heat conduction model including Atangana–Baleanu fractional derivative

Accurately characterizing the thermomechanical parameters of nanoscale systems is essential for understanding their performance and building innovative nanoscale technologies due to their distinct behaviours. Fractional thermal transport models are commonly utilized to correctly depict the heat transfer that occurs in these nanoscale systems. The current study presents a novel mathematical thermoelastic model that incorporates a new fractional differential constitutive equation for heat conduction. This heat equation is useful for understanding the effects of thermal memory. An application of a fractional-time Atangana–Baleanu (AB) derivative with a local and non-singular kernel was utilized in the process of developing the mathematical model that was suggested. To deal with effects that depend on size, nonlocal constitutive relations are introduced. Furthermore, in order to take into consideration, the viscoelastic behaviour of the material at the nanoscale, the fractional Kelvin–Voigt model is utilized. The proposed model is highly effective in properly depicting the unusual thermal conductivity phenomena often found in nanoscale devices. The study also considered the mechanical deformation, temperature variations, and viscoelastic characteristics of the functionally graded (FG) nanostructured beams. The consideration was made that the material characteristics exhibit heterogeneity and continuous variation across the thickness of the beam as the nanobeam transitions from a ceramic composition in the lower region to a metallic composition in the upper region. The complicated thermomechanical features of simply supported viscoelastic nanobeams that were exposed to harmonic heat flow were determined by the application of the model that was constructed. Heterogeneity, nonlocality, and fractional operators are some of the important variables that contribute to its success, and this article provides a full study and illustration of the significance of these characteristics. The results that were obtained have the potential to play a significant role in pushing forward the design and development of tools, materials, and nanostructures that have viscoelastic mechanical characteristics and graded functions.

devices, and various other applications reliant on battery power.The advantageous characteristic of low power consumption is particularly significant in applications like remote sensing and monitoring, where there may be limitations on available power resources 23 .Although nanoscale NEMS devices provide some benefits, they also pose obstacles concerning their manufacture, integration, and dependability.Researchers and engineers persist in exploring novel manufacturing processes, materials, and design approaches to tackle these problems and unleash the full capabilities of nanoscale NEMS technology.These endeavors have resulted in notable progress in developing durable and dependable nanoscale NEMS devices for various applications, such as bioelectronics, nanorobotics, and quantum computing, among others.
Nonlocal elasticity theories have a crucial impact on the examination and development of sophisticated materials, structures, and devices, particularly when dealing with tiny length scales where conventional continuum mechanics theories may not be suitable.These theoretical frameworks offer a more profound comprehension of the complicated mechanical phenomena exhibited by materials possessing intricate microstructures and interfaces 24 .This, in turn, aids in advancing novel materials and structures that exhibit improved performance and dependability.
Two significant factors that alter the mechanical characteristics of nanobeams are the small-scale effect and the surface energy effect 25 .Previous research has investigated these phenomena through the utilization of experiments conducted at the nanoscale and in atomistic models.Moreover, the investigations above have provided evidence indicating that classical continuum mechanics models cannot capture any of these impacts.Consequently, researchers have proposed numerous non-classical continuum models over time to effectively represent the size-dependent characteristics of small structures.These models include the strain gradient Mindlin's elasticity theory 26 , the strain-driven Eringen's nonlocal model 27,28 , the stress-driven Romano and Barretta's model 29 , and other coupled concepts such as the nonlocal strain gradient theory constructed by Lim 30 .Barretta and Marotti de Sciarra 29 also introduced the two-phase local/nonlocal stress and strain gradient model.Romano et al. 31 have extensively debated the potential for mathematical errors in Eringen's model.The primary cause of this is the inherent mismatch between the equilibrium requirements and the higher-order constitutive boundary conditions.
Thermoelasticity is a viable approach for analyzing materials and structures subjected to modest temperature gradients, and it finds extensive use across many engineering disciplines.Usually, thermoelasticity is based on the Fourier heat conduction equation and the idea that temperatures should be equal in a certain area 32 .This theory explains the interaction between temperature gradient and mechanical deformation within materials, taking into account the effect of temperature changes on material hardness and thermal expansion 33 .Generalized thermoelasticity models [34][35][36][37][38][39][40][41][42][43] encompass the fundamental principles of conventional thermoelasticity while incorporating additional considerations, such as the limited wave speeds of thermal disturbances and the non-instantaneous response of materials to temperature variations.These theories look at how the thermal and mechanical fields interact with each other.This helps us get a better picture of the transient thermoelastic phenomena that happen in materials when the temperature changes quickly or when they are loaded and unloaded quickly.
Moore-Gibson-Thompson (MGT) thermoelasticity theory is a new theoretical framework that extends traditional thermoelasticity theory by including the effect of the concept of thermal relaxation times on the transmission of thermal waves within elastic materials.Quintanilla 44,45 formulated the hypothesis above with the aim of offering a more all-encompassing depiction of the transitory thermal and mechanical phenomena seen in materials when exposed to abrupt alterations in temperature.In recent years, the number of studies focusing on developing the theory of thermoelasticity (MGT) has witnessed significant growth [46][47][48][49][50] .Later studies have changed the extended thermoelastic theory into fractional ones by adding different time-fractional derivatives to hyperbolic heat transfer and mass diffusion equations 51,52 .The growing number of fractional calculus applications in both science and engineering served as the impetus for this extension.Due to the vast availability of energy sources that can be readily collected, the micro-and nanoscale vibration-based piezoelectric energy harvester has rapidly become an essential branch of the major emphasis area in modern times.Regarding piezoelectric nanostructures, thermoelastic diffusion, viscoelastic composite structures, and time-fractional order strain, references [53][54][55][56][57][58][59] also provided a comprehensive assessment of current accomplishments and fundamental formulations in the field.
The development of the nonlocal thermo-viscoelastic fractional order model for functionally graded nanoscale beams is a notable progress in the realm of advanced materials and nano-mechanics.This model offers researchers and engineers a potent tool to explore and create functionally graded materials with tailored features.Consequently, it aids in creating cutting-edge and durable nanoscale gadgets and structures for diverse engineering purposes.
The main goal of this work is to study the vibrational behavior of functionally graded (FG) nanoscale viscoelastic beams under a variety of different environmental conditions.We used the theory of non-local elasticity to capture non-local effects on a small scale, and we considered fractional calculus to calculate memory effects.We also used the Kelvin-Voigt fractional viscoelastic model to describe how nanobeams behave in terms of viscoelastic properties.This model is based on the Atangana and Baleanu (AB) differential operators.The Bernoulli-Euler beam theory, a classical beam theory ideal for thin beams that experience modest deformations, was also considered to model the nanobeams investigated in this work.The non-local Bernoulli-Euler beam theory, the fractional Kelvin-Voigt viscoelastic model, and the theory of thermoelasticity are all combined in this study.This gives a complete picture of how FG nanoscale viscoelastic beams behave when they are heated and stressed.Taking into account the effects of small-scale, non-local, and relaxation time, the study explores the effect of different final conditions on the vibrational response of the beam.Finally, the governing equations incorporate a new fractional-order time derivative with a non-singular and non-local kernel.This inclusion allows a more precise depiction of the complex dynamics exhibited by the FG nanopackage as well as its memory-dependent properties.
The proposed model is used to analyze the thermomechanical interactions of nanobeams made of functionally graded (FG) thermo-viscoelastic material that are simply supported and exposed to harmonic thermal flux.The equation for transverse vibration is obtained by the use of Euler-Bernoulli beam theory and Hamilton's principle.A system of governing equations was established and solved utilizing the Laplace transform method in accordance with the given issue.Numerical findings are provided to analyze how the fractional order parameter, nonlocal parameter, and power law index affect the physical variables in the system.The findings from this study may be used to identify and describe different nanostructures, such as nano-electromechanical systems (NEMS), nano-actuators, and others.Studying how functionally graded nano-scale viscoelastic beams vibrate offers important information for designing, optimizing, and analyzing the performance of these nano-structures in many engineering applications.

Formulation Fractional Kelvin-Voigt model
In the context of homogeneous and isotropic substances, the constitutive equations and strain-displacement relations can be formulated as follows 43 : where and µ are the Lame's constants, δ ij is the Kronecker delta, τ ij are stress components, e ij are strain com- ponents, u i are displacement vector's component, γ = (3 + 2µ)α t is the thermal elastic coupling, α t represents the coefficient of linear thermal expansion, θ = T − T 0 is the temperature variation, the variable T represents the absolute temperature, and T 0 denotes the temperature of the medium in its inherent condition.The Young's modulus ( E ) and the Poisson's ratio ( ν ) are connected to the Lamé parameters, and µ , in the following ways: The equation of motion governs the dynamic response of a thermoelastic body to external forces (F i ) and temperature variations.The equation of motion for a thermoelastic body can be expressed under the assumptions of isotropy and homogeneity, taking into account that the deformations are of relatively small magnitude, as follows 43 : The classical Kelvin-Voigt model, a widely recognized viscoelastic model, is employed to characterize the mechanical response of materials that have combined viscous and elastic characteristics.In order to take into account a substance's viscoelastic properties, especially within the Kelvin-Voigt model, the traditional Young's modulus ( E ) is changed to include the effect of viscosity.This model is commonly used to represent the viscoe- lastic behaviour of materials, where both the rate of deformation influences their mechanical properties over time and the duration of applied loads.The modified Young's modulus is expressed as 60 : where E 0 is the elastic Young's modulus and τ v is the vicosity or viscous damping coefficient.Substituting Eqs.(3) and (5) into Eq.(1) leads to 61 : The fractional Kelvin-Voigt model is an expansion of the conventional Kelvin-Voigt model, widely employed for characterizing the viscoelastic properties of materials.A fractional Kelvin-Voigt viscoelastic beam in this study displays viscoelastic properties and is characterised by the fractional Kelvin-Voigt model.Adding fractional calculus to the Kelvin-Voigt model makes it better at showing how viscoelastic materials are, especially when the material behaves in complicated ways, like when it has traits that are passed down through generations or when it remembers things over time.
The constitutive equation governing the behavior of the fractional Kelvin-Voigt model can be expressed as follows 62 : In Eq. ( 7), D α t is fractional differential operator of order α , ( 0 < α < 1 ).Among the many definitions of fractional derivatives, the definition is often used in the Riemann-Liouville sense, which is defined using the relation 7,8 : (1) (5) Vol Recently, other formulations of fractional derivatives involving non-singular kernels have been proposed, such as the Caputo and Fabrizio 10 , Atangana and Baleanu 9 fractional operators.Standard fractional derivatives can be used with these fractional derivative operators, which allows systems that have genetic and memory effects to be represented and studied.The importance of fractional derivatives with non-singular nuclei derives from the observation that some models of dissipative events cannot be adequately described by conventional fractional operators.
In the present work, the Atangana and Baleanu 9 fractional derivative operators will be considered, which is defined as follows: where µ α = α (1−α) and E α is the Mittag-Leffler function.

Nonlocal theory of elasticity
The nonlocal integral theory of elasticity looks at the stress field at a certain point in space by combining the elastic strain field with the right averaging kernel and doing an integral convolution.When applied to fields that change slowly over a characteristic distance, Eringen 28,63 showed that integral-type constitutive equations can be boiled down to partial differential equations.The Eringen model 28 states that the differential constitutive relation can be used to describe the nonlocal constitutive behaviour of a Hookean solid as follows: where σ ij is the nonlocal stress tensor, ∇ 2 is Laplacian operator and ξ = (e 0 a) 2 is the nonlocal parameter in which e 0 represents the constant of the material and e 0 represents the internal characteristic lengths.

Moore-Gibson-Thompson thermoelastic theory
In the context of classical thermoelasticity, the fundamental set of field equations encompasses the classic Fourier law of heat conduction.When these equations are simplified to the displacement-temperature field equations, they exhibit characteristics of a hyperbolic-parabolic nature.This means that a traditional thermoelastic body's response to thermomechanical stress has a propagation velocity that is not limited.Since 1967, a collection of theories known as generalised thermoelasticity theories have been formulated [34][35][36][37][38][39][40] .The models proposed by Green and Naghdi 37,38 have garnered significant interest in the field due to their contributions to generating three distinct forms of generalised thermoelasticity theories.This has subsequently motivated scholars to pursue further investigations in this particular area.Fourier's law has been redefined according to the GN-III model as 37 In this context, the symbol K is used to represent the thermal conductivity of the material, K * is the rate of conductivity, q is used to represent the heat flux vector and V symbolizes the thermal displacement in which V = θ.
The GN-III model similarly predicts the unbounded propagation of thermal waves, consistent with Fourier's conventional law 64 .Hence, the introduction of a modified version of Eq. ( 11) has been proposed in order to address this seeming challenge.In a recent study, Quintanilla 44,45 introduced a novel theory of thermoelasticity in which the Moore-Gibson-Thompson (MGT) equation serves to describe the heat conduction equation.The MGT thermoelasticity concept is a novel theoretical framework that extends the concepts of the LS theory and the GN-III model in the field of thermoelasticity.The incorporation of a relaxation factor in Eq. ( 11) has been implemented in this revision, resulting in 44 The symbol τ 0 represents the parameter associated with time relaxation.Numerous studies (see for example 65 ) have obtained experimental values for τ 0 .
The entropy heat equation is a mathematical formulation that characterises the transfer of thermal energy inside a system by using the concept of entropy.A basic understanding of the relationship between entropy and heat transport is provided by this idea, which is often used to look at irreversible thermodynamic phenomena.The entropy thermal equation is constructed from the second rule of thermodynamics, which postulates that the entropy of an isolated system remains constant or increases with time.The entropy-heat (energy) equation can be expressed as follows: By combining Eqs. ( 12) and ( 13), we may obtain the MGT thermo-viscoelastic heat transfer model with fractional order derivatives as follows Pellicer and Quintanilla 66 have demonstrated the uniqueness and instability of solutions in the concept of thermoelasticity, specifically in relation to the MGT equation.In a recent study conducted by Bazarra et al. 67 , the numerical aspects of the MGT thermoelastic theory were examined.The discrete stability property was also demonstrated in this framework.

Formulation of the problem
A thermoelastic solid nanobeam with functionally graded properties in the Cartesian coordinate system Oxyz will be considered in this study.The x-axis is aligned with the axial direction of the beam, while the y and z axes represent the breadth and thickness, respectively (Fig. 1).At the state of equilibrium, the FG beam remains in a condition of zero strain, stress, and uniform temperature T 0 across its entirety.The analysis focuses on the small flexural deflections exhibited by a nanobeam, which has dimensions of length L (where 0 ≤ x ≤ L ), width b (where −b/2 ≤ x ≤ b/2 ), and thickness h (where −h/2 ≤ x ≤ h/2 .The basic equations governing motion and heat transfer will be considered in the context of generalized (non-Fourier) MGT thermoelasticity.The case will be taken into account in which there are neither body forces, nor external loads ( F i = 0 ) nor even external heat sources ( Q = 0).
There is a smooth change in the material properties only in the z-direction thickness of the nanobeam, which is made of an elastic and isotropic functionally graded material.The FG nanobeam's material properties, such as its Young's modulus ( E ), thermal conductivity ( K ), thermal expansion ( α t ), and mass density ( ρ ), are thought to change continuously along the thickness direction ( z-axis direction).The power function P(z) models this variation by relating the values to the volume fractions of the constituents.
Material properties over thickness show variations depending on the volume fraction, which can be expressed as Here, n p stands for the gradient index, which determines the progressive change in material properties that occurs through the thickness of the nanobeam.The variable P c represents the material property of the pure ceramic, while P m represents the material property of the pure metal.It is important to acknowledge that the material properties of the beam under consideration exhibit a metal-rich composition at the lower surface, located at z = h/2 , and a ceramic-rich composition at the upper surface, located at z = −h/2 .When the value of n p is set as zero, the beam will exhibit homogeneity, indicating that the material composition will consist only of ceramic.
The equations of motion are determined using the principles of Euler-Bernoulli beam theory.Based on this theory, the displacement field at any position along the beam can be expressed as follows: where w is the deflection.
In the context of fractional modeling and viscoelastic functionally graded materials (FGMs), the Young's modulus E m of the material can be represented mathematically as follows: (15) P(z) = P m e n p (2z−h)/h , n p = ln P m /P c .In this representation, the utilisation of fractional order modelling makes it possible to provide a more precise description of the behaviour of the viscoelastic material.Additionally, it provides us with a means to examine and speculate on the behaviour of FGM viscoelastic beams when subjected to many distinct stresses.Using Eqs. ( 15)-( 17), the MGT heat transfer Eq. ( 14) for the nanobeam can be derived as follows where In accordance with Eq. ( 18), it is important to take notice that the coefficients n K , n γ , n K * , and n ρC E are estab- lished by considering the characteristics of ceramic and metal materials.
Assuming that there is no transfer of thermal energy between the upper and lower edges of the nanobeam under study, and that its upper and lower edges are effectively thermally insulated.It has been suggested that there is no temperature difference between the top and bottom surfaces of the nanobeam, and as a result, ∂θ ∂z = 0 at z = ±h/2 .It will also be taken into account that the temperature changes in a sinusoidal form along the thick- ness direction ( z direction) of the nanobeam, and this means that the temperature distribution as a function of z can be expressed in the form of a sine function as follows: By substituting Eq. ( 20) into Eq.( 18) and performing integration throughout the nanobeam thickness from −h/2 to h/2 with respect to z , the following equation can be obtained where Equations ( 15) and ( 16) can be used to obtain the foundational equation that includes the dimensionless nonlocal axial stress σ xx , which appears in Eq. (10), based on the nonlocal theory of thermoviscoelasticity as follows: where The equation for the bending moment, denoted as M , can be derived using Eq. ( 23) as in which The application of Hamilton's principle to micro/nano beams results in the derivation of the Euler-Bernoulli beam transverse motion equation, which is utilized to characterize the lateral displacement of the beam.Based on the expanded Hamiltonian principle, the following equation of motion can be obtained: Vol:.( 1234567890) www.nature.com/scientificreports/When Eq. ( 24) is combined into Eq.( 26), the beam motion equation results as follows: The bending moment given in Eq. ( 24) can be represented as follows when applying expression (26): In order to facilitate analysis and calculation, the following variables are presented dimensionless: After using the nondimensional variables (29) in Eqs. ( 21), (27), and ( 28) and removing the prime, it is possible to obtain: where Equations ( 30) and (31) show the system of equations that governs motion and heat transfer and describe transverse oscillations in functionally graded thermoelastic nanobeams.By solving linear partial differential equations, the deviation w and the function , and thus the rest of the areas of the system, can be determined.For this, the Laplace transform method will be used.

Laplace transform technique
The use of the Laplace transform has proven to be a very effective method for solving linear differential equations with constant coefficients.When using the Laplace transform technique to solve initial value problems, it is usual to use initial conditions, which can be imposed as follows: By utilising the Laplace transform method on the governing Eqs. ( 30)- (32), we are able to derive the following results: where By utilising Eqs.(35) and (36) one can derive the following differential equations that describe the functions and w where ( 27) www.nature.com/scientificreports/Equation ( 33) can be rewritten as The variables h 2 n , n = 1, 2, 3 in Eq. ( 41) represent the roots of the equation The solution to Eq. ( 41) can be written generally as follows: where . The transformed solution in the Laplace field of the displacement u can be derived using Eqs.( 16) and (43)  as follows: In order to determine the bending moment M , the values of the functions w and are entered into Eq.( 37) in the following manner: where

Application
When trying to precisely describe the behavior of the nanobeam and forecast how it will react to various loading scenarios, having a solid understanding of the boundary conditions is necessary.These boundary conditions have a substantial impact on the way the nanobeam behaves; specifically, they affect the nanobeam's deformation and the distribution of stress when it is subjected to different loading situations, such as point loads, dispersed loads, or moments.
In the case where both ends of the nanobeam are simply supported, this means that the beam is held in place in a way that allows it to rotate freely and prevents it from moving in any direction perpendicular to its longitudinal axis at the support points.The simply supported boundary conditions for the nanobeam can be expressed as: The heat flux q(x, t) incorporates thermal effects into the setting of nanobeams, which can affect the nanobe- am's temperature distribution, thermal deformation, and structural integrity.In the current work, it is postulated that the initial end of the nanobeam is subjected to a heat flux q(t) that varies harmonically with time.Therefore, the following thermal conditions are taken into account The symbol ω in Eq. ( 49) represents the periodical frequency of the flow of heat.The frequency parameter ω measures the rate at which the heat flow varies with time, providing information on the frequency at which the heat input pattern repeats during a specified period.In a scenario where there is a constant heat flow, ω = 0 will be set.
When the thermal insulation is applied to the second end of the nanobeam located at position x = L , it sig- nifies the absence of heat transmission at that particular end.This particular boundary condition indicates that the temperature gradient at the terminus of the nanobeam is zero, hence leading to the absence of heat transfer across the boundary.In this case, the following thermal conditions can be taken into account: www.nature.com/scientificreports/When applying the Laplace transform to the boundary conditions ( 47)-( 49), the following equations can be obtained After inserting Eq. ( 43) into the previous boundary conditions, the following system of linear equations can be derived: The indefinite coefficients C n ( n = 1, 2, . . ., 6 ) can be found by solving the above system of linear equations.The final step of the solution is to find the inverse Laplace transform.
When applying the Laplace transform method to solve differential equations, the reverse Laplace transform is a critical but challenging phase.Consequently, algorithms for the numerical inverse Laplace transform are frequently employed to compute the outcomes.There have been a great number of numerical inverse Laplace transform methods developed in order to overcome the challenges associated with Laplace transform inversion [68][69][70][71][72] .All the algorithms utilised have produced satisfactory outcomes.Choosing the appropriate algorithm is contingent upon the specific task at hand.Given that each approach has its vulnerabilities, it is advisable to concurrently implement many algorithms in a programme and evaluate the outcomes for comparison.
This work uses the Riemann sum approximation approach together with computer-generated findings to determine the values of altered fields inside the space-time field, employing the Honig and Hirdes algorithm 72 .The Honig and Hirdes methods 72 are a reliable and quick way to find the numerical inverse Laplace transform, which can be used in the area of fractional-order thermoelasticity.The technique has the potential to be an effective tool for solving differential equations of fractional order.

Numerical results
In this part, the suggested model will be used to examine how changes in model materials and geometrical parameters affect the thermomechanical vibration properties of FGM beams, considering how these properties change with size.Table 1 presents the mechanical properties of the composite material comprising ceramic and metal components.The nanobeam is constructed exclusively with ceramic material, specifically alumina, on its upper side, located at z = −h/2 .Conversely, its bottom surface at z = −h/2 comprises a metallic material, specifically aluminum.
The ratios L/h = 20 and b/h = 0.5 will be chosen based on the assumption that the length of the nanobeam is L = 1 × 10 −9 metres for the purposes of numerical calculations.In addition, it will be taken into account that Table1.Values of the mechanical and thermal parameters of the nanobeam.www.nature.com/scientificreports/ the unit of instantaneous time t , as well as the phase delay duration referred to as τ 0 , is picoseconds.It has been shown that the value of the nonlocal index ξ of carbon nanotubes typically ranges from 0 to 2.0 nm 73 .Numerical calculations of the non-dimensional physical field variables were performed in the case where the length of the nanobeam changes while the variables z = h/6 and t = 0.12 remain constant.The graphical representations of the numerical computations for the bending moment M , temperature change θ , axial displacement u , and lateral vibration w are depicted in Figs. 2, 3, 4 and 5.These figures aim to investigate the impact of the viscosity τ v , non- local parameter ξ and the fractional-order differential actuators α on the characteristics of these fields.The analy- sis of numerical computations and graphical representations can be categorized into three distinct case studies.

Case I: effect of differential operators of fractional order
In order to represent real-world systems characterised by memory effects, the Liouville-Riemann (RL) and Atangana-Baleanu (AB) fractional derivatives offer unique perspectives on fractional derivatives and have distinct advantages and practical applications.Determining which of these alternatives to implement often depends on the specific requirements of the given issue and the formulation's compatibility with the environment of the application.The AB fractional derivative is a variant of the Caputo fractional derivative that integrates the Mittag-Leffler function to model memory impacts accurately.The AB fractional derivative has garnered considerable attention owing to its enhanced capacity to represent memory effects inside certain physical systems correctly.This section will examine the influence of the AB fractional derivative on the responses of the fields under investigation.The numerical findings and comparisons between variables RL and AB will be visually depicted in graphs 2-5.It is worth noting here that in the case of classical viscoelastic theory with integer derivatives (the conventional Kleven-Voigt model), it is possible to set α = 1 .Otherwise, in the case of modified viscoelastic theory with fractional derivatives (fractional Kleven-Voigt model), α = 0.85 and α = 0.75 will be taken into account.In this particular instance, the varying frequency coefficient of the applied heat flow remains fixed at ω = 5 .In addition, we take z = h/4 , ξ = 03 , and τ 0 = 0.02 .The results of this study show that the AB fractional derivative has a big effect on the responses of many domains, especially those with complicated systems that  show memory effects.Also, it is important to point out that the Atangana-Baleanu operator has proven to be an indispensable tool for resolving various dynamical systems, including those that the RL operator cannot manage.
Figure 2 illustrates the impact that fractional-order derivatives have on the dimensionless lateral vibration ( w ) behaviour of viscoelastic simply supported nanobeams in the direction of the axial distance x .The transverse deflection, w , is examined in relation to the distance x while considering various fractional differential operators.Figure 3 shows the effect of fractional actuators on the temperature change θ .Also displayed were the bending moment M and displacement curves u with fractional differential actuators (RL and AB).As shown in Figs. 3  and 4, respectively, it is evident from these figures that modifying the fractional operator, as well as the fractional order parameter, has considerable effects on the behaviour and magnitudes of the various distributions.These changes can be seen in all areas studied.In the same vein, it was discovered that the fractional order coefficient α affects the absolute values of the analysed physical variables, such that a higher value for the fractional order indicator results in higher absolute values.The viscous properties of the nanobeam material may be the main reason for this.It is important to take note of the fact that the deflection w vanishes at the limits of the nanobeam in every scenario; this indicates that it satisfies the boundary requirements that were set for the problem.
When the fractional parameter α changes from α = 0.75 to α = 1 , it is obvious that the value of the solu- tions as determined by the fractional derivative operator AB increases.In addition, the results showed that the fractional factor AB has a reducing effect on the magnitudes of field variables and the propagation of mechanical and thermal waves.The figures show that there is a high degree of correspondence between the solutions in the two cases when the values of the order of the fractional derivatives α are close to 1.This indicates that the results obtained from the two forms of rational operators (RL and AB) are consistent when the fractional order approaches the correct order 61,62 .
In addition, when the results of this study are compared to those obtained using the fractional RL derivative, it can be seen that the solutions change very gradually when they are exposed to the fractional AB derivatives.The selection of an appropriate fractional derivative and fractional order is of utmost importance in the modeling process since it can significantly impact the outcomes.Furthermore, it is evident that the procedure for computing  www.nature.com/scientificreports/ the differential equation using the AB derivative is very straightforward and advantageous, a characteristic that is absent in other forms of fractional derivatives 51,52 .The results of this study can be used to learn more about the physical properties of fractional viscoelastic thermoelastic models with fractional differential AB operator.

Case II: sensitivity to non-local properties
When designing and analysing structures, it is essential to recognise and take into consideration the sensitivity to nonlocal characteristics.Engineers investigate how materials and structures behave under various circumstances and determine how variations in nonlocal characteristics can affect the overall performance using theoretical models, simulations, and experiments.In the field of structural engineering and mechanical sciences, nonlocal factors frequently denote characteristics or circumstances that exert an influence on the performance of a system across a continuum rather than at a discrete location.The consideration of nonlocal characteristics can exert a substantial influence on the general efficiency and reactivity of a structure.One of the goals of this section is to look into what happens to nano-elastic materials when the non-local theoretical flexible modulus ( ξ ) changes.The rationale behind this phenomenon is the significant influence that this parameter exerts on the mechanical and thermal characteristics of flexible nanobodies.The comprehension and interpretation of these non-local impacts are of the utmost importance in order to accurately describe and construct nanostructures.This is because standard continuum mechanics may not possess the capability to properly comprehend the intricate behaviour observed at these minuscule scales.In order to maintain the dependability and efficiency of elastic nanomaterials and devices, it is imperative for engineers and scientists in the field of nanotechnology to consider and include non-local influences 23,24 .
The influence of the non-local component ξ on the dynamic behaviour of the viscous nanobeams is depicted in Figs. 6, 7, 8 and 9.It is important to highlight that in the given scenario, where ξ = 0 , it corresponds to the conventional case, specifically the local viscoelastic model.In contrast, when considering the values ξ = 0.001 , ξ = 0.002 , ξ = 0.003 , and ξ = 0.004 , it indicates the inclusion of the fractional order non-local viscoelastic theory.In the numerical computations conducted for this particular case study, the periodic frequency parameter of the applied heat flow is held constant ( ω = 3 ), as well as other relevant constants such as α = 0.75 , τ 0 = 0.02 ,   29,74 .
The results in Fig. 6 show that the non-local elastic parameter ξ plays a large role in knowing how the deflec- tion w is distributed in nanoelastic structures.Non-local effects therefore become important at the nanoscale, where the dimensions of the material are very small.This makes the usual assumptions of classical elasticity less valid.Changes in the non-local parameter ξ could cause changes in the way deflection w is distributed, which could affect the stability and integrity of nano-elastic structures.
The data in Fig. 7 shows that non-local processes (non-local elastic parameter ξ ) can change the way tem- perature change ( θ ) in nanoelastic structures.When looking at thermal properties at the nanoscale, there is a lot of complexity.For example, the non-local elastic factor ξ can have a big effect on how heat moves through or leaves a material.Changes in the non-local parameter ξ could cause changes in the temperature profiles, which would have an effect on the thermodynamic properties and performance of nano-elastic systems 18,20 .
The analysis of Fig. 8 shows that the elastic beams and nanostructures can deform in ways different from the conventional case, where the non-local elastic parameter affects the size and shape of the deformation ( u ).Any change or addition to the non-local parameter ξ can cause changes and improvements in the amount of defor- mation u exhibited by the material, which in turn affects its mechanical response and structural properties.This finding holds great importance since understanding deformation and its regulation at the nanoscale are crucial in the process of creating and reinforcing nanostructures designed for practical purposes.It is also clear from the curves in Fig. 9 that the effect of the non-local theoretical elastic modulus on the bending moment ( M ) variation pattern is essential in the context of elastic nanomaterials and cannot be neglected.
In certain instances, the consideration of nonlocal characteristics can be effectively managed by employing optimisation methodologies, which aim to determine the optimal values of parameters that yield the desired  www.nature.com/scientificreports/structural performance.This process may entail modifying material qualities, geometric configurations, or other variables in order to attain the necessary equilibrium between strength, stability, and additional performance objectives 22,23 .

Case III: the effect of viscosity parameter
Nanobeams, characterised by their nanoscale dimensions, possess distinct mechanical properties that deviate from those observed in macroscopic structures.The addition of viscous damping ( τ v ) is essential in appropriately describing the dynamic behaviour of nanobeams as they undergo vibrational motion.There is a big effect of the viscous damping coefficient τ v , which is also known as the viscosity parameter, on how elastic nanobeams vibrate.Viscous damping ( τ v ) is a mechanical phenomenon that converts mechanical energy into heat, resulting in significant impacts on the dynamic behaviour of objects at the nanoscale.The addition of viscous damping changes the resonance characteristics of nanobeams.The stiffness properties of the material and the characteristics of damping are two important factors that affect the phenomenon of resonance.A thorough examination of damping is important in order to accurately forecast and regulate resonance events in nanobeam constructions.It will be looked into how the viscous damping factor ( τ v ) affects the vibrational properties of the Euler-Bernoulli flexible nanobeam in the last scenario of testing and evaluation.
It is noteworthy that in the case of an elastic material, where the viscosity parameter ( τ v ) is equal to zero, the outcomes of the elastic nonlocality concept are demonstrated.In the numerical calculations in this case, it was taken into account that the other effective constants are constant ( α = 0.75 , ξ = 0.003 , τ 0 = 0.02 , and ω = 3 ).loss of vibrational energy, which in turn leads to a greater amplitude of deflection and distortion 47,50 .This is of utmost significance in nanomechanical systems, as accurate control of deflection is essential for their operation 56 .
Figure 13 clearly shows how the viscous damping ( τ v ) changes the way bending moments are distributed in the nanobeam when it is dynamically loaded with harmonic heat flux.As the level of damping increases, both the amplitude and duration of the bending moment can decrease.The comprehension of the structural integrity of nanoscale systems under dynamic forces is of utmost importance.
The dynamic stability of the nanobeam is significantly influenced by the viscosity and damping parameter τ v .The relationship between stiffness, damping, and viscosity impacts whether the amount of vibration remains limited or expands exponentially over time.Comprehending its inherent stability is crucial in order to forecast the extended-term performance of the nanobeam.In conclusion, a crucial factor in understanding how elastic nanobeams vibrate is the viscous damping coefficient 59,60 .For the purpose of designing and optimising nanomechanical systems and devices, it is imperative that their incorporation into mathematical models be used in order to precisely forecast and comprehend the changing behaviour of viscoelastic nanostructures 61 .

Conclusion
The current study presents the incorporation of the Kelvin-Voigt fractional viscoelastic model into the FG Euler-Bernoulli nanobeam, utilising the nonlocal elastic theory and the Atangana-Baleanu (AB) fractional derivative.This approach aims to characterise the dynamic response of nanostructures, encompassing viscoelastic properties and the influence of small-scale phenomena.We use this model to study and develop functionally graded nanomechanical systems because it is crucial to accurately predict their viscoelastic response.We assume that the characteristics of the nanobeam gradually vary across its thickness, transitioning from a purely ceramic composition to a purely metallic composition.This study examines and analyses the impact of many factors, including fractional order, damping parameters, and nonlocality, on the vibrational characteristics of nanobeams.Based on the results of this research, the following conclusions can be emphasised:  www.nature.com/scientificreports/ • By adding fractional derivatives to the Kelvin-Voigt model in the context of fractional derivative viscoelastic- ity, we can represent the memory and heredity effects in the material's response.• The Kelvin-Voigt fractional derivative viscoelastic model describes how the nanobeam moves and reacts to outside forces or moments.It does this by showing how the beam bends and deforms.• The results showed that the fractional factor AB has a reducing effect on the magnitudes of field variables and the propagation of mechanical and thermal waves.• Changes in the non-local parameter could alter the distribution of deflection, which could impact the stability and integrity of nano-elastic structures.• Any change or increase in non-local parameter values can lead to changes and enhancements in the mate- rial's deformation, subsequently influencing its mechanical response and structural properties.This finding holds great importance since understanding deformation and its regulation at the nanoscale are crucial in the process of creating and reinforcing nanostructures designed for specific purposes.• Higher viscosity materials are likely to experience more dramatic temperature changes when viscous nano- beams undergo thermal loading.The change in temperature also has a direct effect on the damping process, which in turn can affect the overall thermal stability of the nanobeam.
Finally, critical factors in understanding how viscoelastic nanobeams vibrate are the viscous damping coefficient and nonlocality.To design and improve nanomechanical systems and devices, we must use mathematical models to accurately predict and understand how viscoelastic nanostructures change over time.

Figure 2 .
Figure 2. The transverse displacement w against x for different fractional differential operators.

Figure 3 .
Figure 3.The temperature change θ against x for different fractional differential operators.

Figure 4 .
Figure 4.The axial displacement u against x for different fractional differential operators.

Figure 5 .
Figure 5.The bending moment M against x for different fractional differential operators.

Figure 6 .
Figure 6.Variation of tangential displacement w versus non-local parameter ξ.

Figure 8 .
Figure 8. Variation of the axial displacement u versus non-local parameter ξ.

Figure 9 .
Figure 9. Variation of the bending moment M versus non-local parameter ξ.

Figure 12 .
Figure 12.The axial displacement u versus viscous damping coefficient τ v .

Figure 13 .
Figure 13.The bending moment M versus viscous damping coefficient τ d .